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ABSTRACT 

We demonstrate that the process e^e~ tth at the NLC 
provides a powerful tool for extracting the tt (Yukawa) cou- 



have 



plings of the h. In combination with the e 



Zh pro- 



cess, an accurate determination of the ZZ coupling of the 
h is also possible. The resulting ability to distinguish dif- 
ferent models of the Higgs sector is illustrated by detailed 
studies for two-Higgs-doublet models. 

In extensions of the Standard Model (SM) there are mul- 
tiple neutral Higgs bosons. Their masses and couplings 
are often dependent upon many parameters; CP- violating 
mixing of CP-even with CP-odd neutral Higgs fields is gen- 
erally possible. Thus, if Higgs boson(s) exist and are dis- 
covered at future colliders, it will be extremely important 
to determine both the magnitude and the CP nature of 
their couphngs[0 ||, H). 

It has been shown that the process e+e^ tth [h is our 
notation for a generic neutral Higgs boson) at the proposed 
Next Linear Collider (NLC) can provide rather accurate 
determinations of the CP-even and CP-odd tth Yukawa 
couplings and at least a rough value for the ZZh coupling 
1^. The e+e" Zh process provides a direct measure- 
ment of the Z Zh coupling when analyzed in the missing- 
mass mode with Z — > e"'"e~, ^'^ fi~ . Here, we demonstrate 
that the accuracy with which the couplings can be deter- 
mined using these two processes will discriminate in a very 
decisive manner between different models for the Higgs bo- 
son sector. By way of illustration, we will consider the 
general Two-Higgs-Doublet Model (2HDM), for which the 
masses and couplings of the three neutral Higgs bosons 
are all free parameters. For simplicity, we focus on the h 
that is the lightest of the three eigenstates, and assume 
that it is sufficiently light compared to the other two that 
e+e~ tth is not sensitive to the other Higgs states. 

The relevant Feynman rules for the Zh and e"'"e~ — > tth 
processes can be parameterized as: 



tth : —t{a + ibjQ)t 



gmt 
2mw 



ZZh 



gmz 
cos{9w 



where g is the usual electroweak coupling constant, 
the SM, 

a= 1 , 6 = 0, c= 1 . 



(1) 



For 



(2) 



For the 2HDM, the couplings are more complicated. We 
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a ~ . , b — Rsj cot /? , c — Rij cos /3 + i?2j sin /^5) 

where j = 1,2,3 indicates one of the three Higgs mass 
eigenstates, tan/3 is the ratio of the vacua of the neutral 
members of the two Higgs doublets (we assume a type-II 
2HDM), and Rij is a 3 x 3 orthogonal matrix which spec- 
ifies the transformation between the 2HDM Higgs fields 
and the Higgs boson mass eigenstates. We employ the pa- 
rameterization: 



R = 




S1C3 

C1C2C3 - S2S3 
-C1S2C3 - C2S3 



S1S3 

C1C2S3 + S2C3 
-C1S2S3 + C2C3 



(4) 



where Si — sina^ and Ci — cos a,;. Without loss of gener- 
ality, we identify the lightest Higgs h with the j = 1 mass 
eigenstate. In this case, we have 



5lC2 

sin f3 



b — siS2Cotp, c = ci cos/3 — S1C2 sin/3 



(5) 

The h has CP-violating couplings if either ah or be 0. 

We make a few remarks based on Eq. (||) regarding spe- 
cial limiting cases. 

• In the 2HDM context one can always reproduce the 
SM couplings of Eq. (||) for any given tan (3 by taking 
ai = P and a2 = tt. If tan/3 is not large then deter- 
mination of tan /3 would be possible via observation of 
one of the other Higgs bosons {j = 2,3). However, if 
tan /? is large then ai 7r/2; coupled with a2 ''^ tt, 
this implies that the remaining Higgs bosons (j = 2, 3) 
will have small tth couplings, cx 1/ tan/3, and would 
not be easily probed via the tth final state. 

• At large tan/3, 6 '--^ and a ~ c ~ — S1C2 (unless 
S1C2 as in the SM limit); sensitivity to the ex- 
act value of tan /3 is lost and ai and a2 cannot be 
independently determined. 

• For very large tan/3 ^ mt/mb, the bbh couplings can 
be large, in which case the e^e'' — > bbh rate would be 
significant and could be analyzed using the procedures 
to be discussed here for the tth final state. 

In this report, we will focus on models with tan/? in the 
vicinity of 1 in order to display the full potential of the tth 
final state. 

The e^e~ — > Zh differential cross section is kinemati- 
cally trivial. The only useful observable is the total cross 



section, ariZh), which is proportional to c^. The differ- 
ential e^e~ — > tth cross section, without measuring the 
polarizations of the fermions, contains five distinct terms: 



ci — ; C2 — ; C3 ~ be; C4 — ; — ac, 



(6) 



and the fi{(t>) are theoretically known functions of the 
Higgs mass , the machine energy ^/s, and the final state 
phase space variables, 4>, but are not dependent on the 
model. Note the absence of any term proportional to ab. 
The C3 = be term is explicitly CP-violating. The total 
cross section, aTitih), is a particular linear superposition 
of the a: 



aritth) 



da 



-d(j) = cigi + C2.g2 + C4,g4 + csffs 



(7) 



where the gi are functions of rrih and y/s, and specific ex- 
perimental cuts. There is no contribution from the CP- 
violating C3 = be component of da/d(f>. Since many differ- 
ent Higgs parameter choices can yield any given value of 
<tt {tth) , it is vital to make use of the much greater infor- 
mation embodied in the detailed dependence of da /d(j) on 
the (j) variables. 

The statistically optimal technique for extracting the 
Ci using da/d(j) was developed in Ref. One employs 
weighting functions Wi{4>) such that / Wi(4>)[da / dcj)] = Ci, 
where the Wi[4)) are uniquely defined by demanding that 
the statistical error in the determination of the Ci is min- 
imized; this is the choice such that the entire covari- 
ance matrix is at a stationary point with respect to vary- 
ing the functional forms for the Wi{(j)) while maintaining 
/ Wi[(p)fj{(j))d(j) — Sij. The weighting functions are given 
in Ref. f^. By employing them, one finds 

Ci = ^ Xikh = ^ ^^ifc^4 , where 4 = / fk{<t')d(l) , 

k k 

(8) 
(9) 



with 



ik 



[da/d(t)] 

If there are experimental cuts that exclude a portion of 
the phase space in (jj, they should be included in computing 
Mik via the / d(j) appearing in Eq. (||) in order that opti- 
mal statistics be achieved in the presence of the cuts. Since 
the cuts that will be employed are detector dependent and 
cannot be determined at this time, we have opted to com- 
pute Mik in the examples to follow without including any 
cuts. However, we will reduce the total event rate by an 
overall efficiency factor, the magnitude for which will be 
chosen so as to reflect a reduction due to cuts. 

The covariance matrix corresponding to Mik of Eq. 

is 



Vij = {AciAcj 



Mr^^ aritth) 



(10) 



N{tth) ' 

where N(tth) — Lcs{tth)aT(tth) is the total number 
of events, with Lcsitth) being the effective luminosity: 



Lcs{tth) — e{tth)Ltotsd, where itotai is the total inte- 
grated luminosity and e(tth) is the efficiency, including 
branching ratios for the tth to decay into the useful fi- 
nal states. Since identification of the t and t requires 
that one decay semi-leptonically and the other hadroni- 
cally, e{tth) < 2B{t lvb)B{t 2jb) - 0.44. Depending 
upon how the h decays, there may be a further loss for fo- 
cusing on reconstructable h final state decays. There will 
also be cuts and detector efficiencies. We adopt the value 
of e{tth) — 0.1. For the m^j = 100 GeV value that we shall 
focus on, for which B{h bb) ~ 0.9 is likely, this is fairly 
conservative. 

In order to compute the expected experimental errors for 
the Ci, we first compute Alik [using Monte Carlo integra- 
tion in Eq. without cuts] and thence, via Eq. (10), the 
covariance matrix V for the given input model. The confi- 
dence level with which one can rule out parameter choices 
different from those of the input model is then determined 
by the associated value: 



X^tth) 



E 



(c,-cO)(c,-c°)V^- 



with 



1 ^ M^.Nitth) 

ij 



(11) 



(12) 



aritth) 

In Eq. (|l^), the c° are the values for the input model and 
the Ci are functions of the location in a, b, c parameter space 
of the alternative model, see Eq. (||). Sensitivity of x^{tth) 
to the a, b, c parameters is thus directly determined by the 
covariance matrix for a given model. Typically, one finds 
that sensitivity to ci — is largest, while the weakest 
sensitivity is to C3 = be. 

We note that x^itih) implicitly includes a contribution 
due to the difference in aritth) for the input model as com- 
pared to the alternative models. In what follows, we shall 
be implicitly assuming that the only errors in aritth) are 
the statistical ones as incorporated in x^{tth) in Eq. (p]). 
However, we note that aritth) will be subject to system- 
atic error as well. The main uncertainties arise from the 
fact that one must observe iih production in one or more 
particular final states, leading to uncertainty in aritth) to 
the extent that the t and h branching ratios and/or the 
detection efficiencies for these particular final state(s) are 
uncertain. Thus, we will be implicitly assuming that these 
uncertainties can be kept below the level of the simple sta- 
tistical uncertainty. 

The statistical analysis for the e+e~ Zh process 
is completely straightforward. A direct (i.e. indepen- 
dent of Higgs branching ratios) measurement of C4 — (? 
is obtained when the h is isolated via a peak in the 
[(Pe+ + Pe- — Pz)^]^^^ missing mass distribution, where 
we require Z e+e~, /i^/i" in order to be assured of 
the cleanest possible analysis and most reliable absolute 
normalization. The number of Zh events is given by 
NiZh) = L,siZh)ar{Zh), where L.siZh) = e(Z/i)Ltotai, 
with e{Zh) being the efficiency for detecting the events 



using Z — !■ e+e^,/i+/i^ decays: e{Zh) = B{Z — > 
e'^e~ , n'^ fi~)i{Zh). We take e{Zh) — 0.5 for the rem- 
nant efficiency associated with cuts and overall detector 
efficiencies. The relative accuracy of the measurement of 
C4 is simply given by \/ ^ N{Zh), and thus the associ- 
ated with choosing a value of C4 that differs from that of 
the input model value is given by 

^^(Zh)^^^±^^N{Zh). (f3) 

The total associated with choosing values for a, b and c 
that differ from the input model values is given by summing 
the tth and Zh results: 

=X^(tth) + x\Zh). (14) 

We now provide several examples. We take ruh = 
100 GeV and = ITeV. We assume Ltotai = 500 fb"^ 
(as achieved for 2 1/2 years of running at Lyoar — 
200 fb"^). This gives Lesltth) = 50 fb"^ and LcB{Zh) = 
16.9 fb^^. Assuming no cuts, the above to/j and ^/s imply 
ariZh) = 13.6cl fb, yielding N{Zh) = 229. Sc^ events for 
a given input model. For crTitth) [see Eq. (Q)] we find (fb 
units) 

51 ^ 2.70 , 32 = 0.530 , 54 = 0.083 , 35 = -0.055 . 

(15) 

Note the insensitivity of axitth) to ac, and very modest 
sensitivity to c^. 
We consider three input model cases: 

• SM: We assume that the input model is such that the 
Higgs has SM couphngs, Eq. (|). From Eqs. (|), (0) 
and (111), we find aritth) = 2.73 fb, yielding N{tth) - 
136 for Leffitih) = 50 fb"\ For LasiZh) = 16.9 fb"^ 
we obtain N{Zh) - 229. 

• 2HDM(I): We assume that the input model is the 
2HDM model with tan^ = 0.5, ai = 7r/4, and 
a2 = 7r/4, yielding a = -1.118, 6 = 1, c = 0.4088. 
In this case, as compared to SM couplings, (Jt{ZK) is 
smaller, yielding N{Zh) ~ 38, and aT{tth) is larger, 
aritth) = 3.94 fb, yielding N{tth) - 197. 

• 2HDM(II): We assume that the input model is the 
2HDM model with tan/3 = 0.5, ai = 7r/4, and 
a2 = 7r/2, yielding a = 0, & = 1.414, c = 0.6325. 
In this case, as compared to SM couplings, <TT{Zh) 
is smaller, yielding N{Zh) ^ 92, and aritth) is also 
smaller, aritth) = 1.09 fb, yielding Nitth) - 55. 

In all cases, we use Monte Carlo integration to compute 
Mik as given in Eq. (||), and Eqs. ( p^ ) and ( |l^ ) to com- 
pute the matrix V~^; all depend upon the input c". In 
computing Mik we do not 

It is useful to note that the la statistical errors (ex- 
pressed in percentage terms) in aritth) and ariZh), cor- 
responding to the above-quoted event rates are: 

SM 2HDM(I) 2HDM(II) 
aritth) : ±8.6% ±7.1% ±13.5% (16) 
ariZh) : ±6.6% ±16.2% ±10.4% 



We believe that the systematic errors in aritth) and 
ariZh) will be smaller than the above numbers given that 
the detector efficiencies and the relevant t and h branching 
ratios should be very well known by the time this analysis 
is performed. This is presumed to be the case in obtaining 
the numerical results that follow. 

The accuracy with which the 2HDM parameters can be 
determined is illustrated in Figs. |l|, ^ and ^ Each figure 
has six windows. In each window of the three figures, a 
filled central region, an empty band, and a filled band may 
all be visible. The central region is the <1 region, the 
empty band is the 1 < ^ 4 region, and the outer filled 
band is the 4 < x^ < 9 region. If no filled central region 
is visible, the central region being empty, then this means 
that x^ < 1 was not possible. If only a completely filled 
region appears, then x^ < A was not possible. In the three 
left-hand windows of each of the three figures, results are 
displayed for the case where the input model is a 2IIDM 
constrained so as to reproduce the SM couplings when 
tan (3 = 0.5, 1.0 or 1.5. In the right-hand windows we show 
results for 2HDM(I) with tan/3 0.5 and tan/3 = 1 and for 
2HDM(II) with tan/? 0.5. For all 2HDM(I) [2HDM(II)] 
parameter choices, x^ > 9 if tan/? = 1.5 [tan/3 = 1.0 or 
1.5]. 

In Fig. |l| we show the above-described regions in 
(ai/TT, Q!2/7r) parameter space. We restrict the plot to < 
ai < TT and < a2 < tt. (Since da/d(f>itth) and ariZh) 
are only sensitive to a^, c^, b^, ac and be, nothing changes 
if we simultaneously flip the signs of a, 6, c. Restricting to 
< ai.2 < TT avoids this ambiguity.) Note the fact that for 
all but 2IIDM(II) the regions in the < ai < 7r/2 domain 
are self-similar to the regions in the 7r/2 < ai < 1 domain 
obtained by ai tt — ai and a2 — 0:2, which changes 
the signs of a and c, but not b. 

From this figure, we observe the following. 

• In the case of SM input couplings, ai_2 must lie 
close to the a2 = tt and ai = /3 values that yield 
a = c= l,6 = 0, or else to the 02 = 0, ai = tt — /3 val- 
ues that yield a = c = — 1,6 = 0, leaving a^, c^, ac un- 
changed. Note that the different x^ regions all shrink 
with increasing tan /3. 

• In the case of 2IIDM(I), only tan/3 = 0.5 (the input 
value) allows x^ < 1, and the x^ < 1 region corre- 
sponds closely to the input ai^2 values of ai = a2 = 
7r/4 or the a,c sign-flipped ai — a2 — 37r/4 values. 
(Self-similarity of the latter regions to the former 
is a consequence of four facts: i) be ^ 0.41 is not 
large; ii) sensitivity of da/dcj) to C3 = be is weak; iii) 

= 1.25 is large; and iv) sensitivity of da/dcj) to 
is substantial.) If we allow 1 < x^ ^ 4, the allowed 
regions expand considerably, and for 4 < x^ < 9 there 
are two more regions that develop with ai_2 values 
that are very different from the input values. Fur- 
ther, we observe that tan/3 = 1.0 would be allowed at 
the x^ > 1 level for yet another region of a 1.2 values. 
Values of tan/3 > 1.5 are excluded at the x^ < 9 level. 
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Figure 1: Regions of < 1, 1 < < 4 and 4 < < 9 in 
the ai/iT (horizontal axis) and 02/7^ (vertical axis) plane. 
See text for details. 



Figure 2; Regions of < 1, 1 < < 4 and 4 < < 9 
in the a (horizontal axis) and c (vertical axis) plane. See 
text for details. 



• In the case of 2HDM(II), only tan/3 = 0.5 (the input 
value) allows x^ < 9. The x^ ^ 1 region corresponds 
closely to the input values of ai — 7r/4 and a2 = n/2. 
An alternative region with ai — > tt — ai develops for 
4 < x^ ^ 9. (Self similarity under ai^2 — *■ tt — 0:1.2 is 



not present since, unlike 2HDM(I), 
0.89 is fairly large.) 



= and be 



Note that in the non-SM 2HDM(I) and 2HDM(II) cases 
we obtain an approximate determination of tan f3. 

The implications for the a, 6, c couplings appear in Figs. ^ 
and|^, in which the x^ regions are plotted in the (a, c) and 
(a, b) planes. All regions with x'^ ^ 9 are shown in the 
(a, c) plane figure. In the case of the (a, 6) plane, only 
a > is shown except in the 2IIDM(II) case where both 
the a > and a < regions are shown (note the difference 
in horizontal axis labelling for this case). For the SM and 
2HDM(I) self-similar region to the one displayed 

for a > is obtained for a < by flipping about the 
= axis. (Regions with 6 < do not emerge.) From the 



figures, we observe the following. 

• For SM input couplings, the output values of a and c 
must be very close to the a = l,c = 1 input values, 
or the alternative a = — l,c=— 1 flip, either of which 
require b > 0. (The x^ ^ 1 regions are very small dots 
in the (a, c) plane; careful examination of the picture 
is required.) The value of b is only moderately well- 
constrained when tan /3 — 0.5, with b < 0.4 (0.7) being 
allowed at the x^ ^ 1 (x^ 1^ 4) level. The constraint 
on b becomes much tighter as tan P increases, with 
b < 0.2 being required for x^ 5; 4 once tan/3 > 1.5. 

• For 2IIDM(I) input, the tan/3 = 0.5 windows of the 
(a, c) and (a, b) planes show that the a, 5, c couplings 
are all very well-determined at the x^ < 1 level (up 
to the sign- flip of a and c). Substantial flexibility in 
b develops for 1 < x^ 5; 4. For 4 < x^ ^ 9, a region 
where a has changed sign (but not c) develops. For 
tan/3 = 1.0, x^ ^ 1 is not possible, but for 1 < x^ !i 4, 
a solution develops that has the wrong sign of ac and 



a very distorted value of b. ^ ^ is not possible for 
tan P = 1.5. 

• For 2HDM(II) input, the a, b, c are again well- 
determined if we demand ^ li ^-nd ^ 4 al- 
lows much less flexibility than in the 2HDM(I) case. 
However, 4 < x^ !i 9 allows a a solution with the 
flipped sign of ac and slightly distorted b values. [In 
the (a, 6) plane 2HDM(II) window, the three different 
X^ regions associated with the correct sign of ac are 
somewhat obscured by the strange extra blob associ- 
ated with 4 < x^ < 9 and the wrong sign of ac] 
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Figure 3: Regions of x^ < 1, 1 < X^ ^ 4 and 4 < x^ < 9 
in the a (horizontal axis) and b (vertical axis) plane. See 
text for details. 

In conclusion, we note that it is very possible (some 
would say probable) that the SM is not correct. In this 
case, and if there is a weakly-coupled Higgs sector, there 
will certainly be Higgs bosons that do not have SM-like 
couplings. This is true even if one neutral Higgs is very 
SM-like (as for example is very probable in the minimal 
supersymmetric model), since the others must have very 



small ZZ coupling and can have all manner of tt cou- 
plings. Thus, it will be crucial to determine if an observed 
Higgs boson fits into a given model context, such as the 
two-Higgs-doublet model, and to determine the model pa- 
rameters and associated couplings for acceptable solutions. 
By doing this for all the Higgs bosons we would be able to 
completely fix the Higgs sector model and parameters. 

In this report, we have examined the possibility of car- 
rying out such a program by applying the optimal analysis 
procedure of Ref. to the e+e^ — > tih differential cross 
section and measuring the e~^e~ —^ Zh total cross section. 
Using Ltotai — 500 fb~^ of data from the NLC operating 
with -y/s = 1 TeV, we have demonstrated that for models 
with a reasonable tih event rate the couplings of a 100 GeV 
2HDM Higgs boson can be determined with substantial ac- 
curacy at the la level. However, for this luminosity some 
ambiguities begin to arise in the 1 — 2cr range. Ambiguities 
at the < 1(7 level could arise if systematic uncertainties in 
the experimental determination of the overall normaliza- 
tion of the tth and Zh total cross sections are not small 
compared to the statistical accuracies. At larger Higgs 
masses, statistics will deteriorate; higher Ltotai will be re- 
quired to avoid significant ambiguity. However, even when 
ambiguities emerge, we have found that they are usually 
sufficiently limited that the type of analysis presented here 
will make a critical contribution to gaining a clear under- 
standing of the exact nature of all the Higgs bosons. Cer- 
tainly, these procedures will provide a powerful means for 
distinguishing between substantially different models. We 
urge our experimental colleagues to carry out fully realistic 
simulations of this type of analysis. 
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